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Photoproduction of neutral-kaon pairs is studied from the perspective of CP and 
CPT studies. Interference of the P and S waves, with the former due to diffrac- 
tive production and the latter to fo/a production, is shown to enhance the CPT 
reach. Results are presented of Monte Carlo studies based on rates expected in future 
experiments. 



1. Introduction. Neutral-meson oscillations provide a sensitive tool for testing CPT 
symmetry M. Impressive bounds have been achieved in experiments with both K 
mesons [0, [!| and mesons Although the CPT theorem |J guarantees that the 
standard model preserves CPT by virtue of its construction as a Lorentz-invariant 
quantum field theory, violations of CPT could be exhibited in a more fundamental 
description of nature that incorporates physics at the Planck scale ||. For example, 
CPT and Lorentz violation may naturally arise in string theory J7| and could lead to 
effects in the various neutral-meson systems ||. It is therefore valuable to identify 
additional possibilities for future experiments with exceptional sensitivity to CPT 
violation. 

In the K system, CPT bounds of a few parts in 10 19 now exist on the ratio of 
the kaon-antikaon mass difference to the kaon mass 0, |3|. Future improvements over 
these bounds are expected from specialized experiments at kaon factories, such as 
KLOE at DAPHNE ||. Indeed, it has long been recognized that high-luminosity 
(f) factories are ideal for studies of CP and CPT violation |Kj because a decaying 
meson produces a well-defined flux of C-odd K S K L pairs. 

In the present work, we study the possibility that interference effects between the 
C-odd wave and a coherent C-even K°K° background could be used to enhance fur- 
ther the sensitivity to parameters describing the weak K° and K° decay amplitudes. 
At an e + e _ collider, the <fi mesons are essentially the only source of K°K° pairs. 
However, such studies potentially suffer from a C-even background from the decay 
lyj 4> —> 7/0 (980) — > / jK°K°. Fortunately, recent measurements by the SND and 
CMD-2 collaborations report fl2j a branching ratio B(0 — > 7/0) = 1-9 — 3.5 x 10 , 
which is insignificant for CP and CPT studies ||13| . 

Here, we consider instead hadronic production, where a different mechanism exists. 
In particular, in photoproduction near the <fi meson peak, a significant S'-wave back- 
ground has been measured and attributed to the decays of the /o(980) and ao(980) 
mesons [14, [L5|]. In what follows, we show how differences in the angular distributions 



of the S and P waves could be exploited for CP and CPT studies. 

Nature is believed to be described at low energies by a quantum field theory. The 
CPT theorem then suggests that any violations of CPT invariance must be accom- 
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panied by Lorentz violation. At the level of the standard model, small violations of 
Lorentz and CPT invariance can be introduced via additional terms in the lagrangian, 



which yields a general standard-model extension Ulb] . In this general context, it turns 
out that neutral-meson oscillations provide a unique sensitivity to a class of parame- 
ters for CPT violation associated with flavor- changing effects in the quark sector [[17] . 



In fact, to leading order in small parameters, tests with neutral mesons are indepen- 
dent of all others performed to date, including comparative measurements in Penning 



traps [181, spectroscopy of hydrogen and antihydrogen | 19fl , measurements of muon 



properties [[2(J, clock-comparison experiments ||21| , observations of the behavior of 



a spin-polarized torsion pendulum [p2 |, measurements of cosmological birefringence 



p3| , [IB| , and observations of the baryon asymmetry f24 |. 

In the K system, experiments have reached a sensitivity of parts in 10 20 to certain 
parameters for CPT violation in the standard- model extension JJ, |3], [L7|]. At leading 
order, four coefficients control CPT violation for kaons. Experiments can therefore 
place four independent CPT bounds, only two of which have been obtained to date. 
The photoproduction of mesons offers a distinct experimental arena with the po- 
tential to bound new combinations of parameters, a possibility which may well merit 
careful investigation. However, with the exception of a few remarks below, the scope 
of the present exploratory work is limited to the demonstration that a background- 
enhanced CPT reach is possible in photoproduction. 

2. Theory for S-P interference. In the reaction jp — > Xp —>■ K°K°p with J PC (X) = 
1 ((p), ++ (/o) or ++ (ao), the K°K° wave function in the rest frame of the pair 
can be written as 

\i) = hj d 2 q [a p (q)\Ks(q)K L (-q))+a s (q) (\K L (q)K L (-q)) - \K s (q)K s (-q))) . 

(1) 

Here, the first term represents the P wave and comes from the decay and other 
coherent odd-parity backgrounds with a p (q) = J2 m a m^im(<?), where q is a unit vec- 
tor. The second term comes from 5-wave and other even-parity backgrounds with 
a s = ^aQY 00 (q). The photoproduction amplitudes a J m = a^(s, t, m K ^) describe the 
dynamics of the production process. Since KK photoproduction is dominated by 
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helicity-nonflip diffraction for which the a l ±1 coefficients are the largest, the subse- 
quent KK evolution is best studied in the s-channel helicity system, with the z-axis 
defined in the direction opposite to the direction of flight of the outgoing nucleon in 
the KK rest frame [p5 1 . In general an incoherent .S-wave background is also present 
but is largely irrelevant for our analysis. We comment on its effect later. 

The normalization h in Eq. (|T]) is chosen so that the KK photoproduction rate is 
given by 

dN KR = F <Hrp-> kk p ) w 

dtdm KR d 2 q dtdm KR 
where F is the photon flux. In this equation, the KK angular distribution W (q) is 
taken to be unit normalized, / dqW(q) = 1, and is given by 

W(q)=J2Wjj,= £ P J S>YjmYj*, m ,. (3) 

JJ' Jm,J'm' 

Here, p^m' = a m a m* are elements of the spin density matrix. Summation over nucleon 
helicities is implicit. 

The two kaons decay into final states f a , a = 1,2, with amplitudes given by 

A{K S {l) -> fa) = a a ,s(L) exp(i(pa,s(L))exp( - im S (L)t a - T a s[L) t a /2), (4) 

where ^s(l) are ^ ne corresponding partial decay widths. For convenience, we define 
as usual for each fixed a the parameter 

r] a = 1 77 a | exp(i(j) a ) = a ajL /a atS (5) 

as the ratio of the amplitude for the transition between f a and Ki to that between 
f a and K s . For the moment, the f a are kept general. 

The production rate R(t\,t2] q) of the final state f a with momentum direction of 
fi specified by the solid angle q is 

R( tl , t 2 , q) = j dm^dt ^^Jl ■ (6) 
This expression can be expanded as 

R(t 1 ,h,q) = N KR (\Mpp\W P p+ £ \At? s \W? s +\M ss \W 3 s), (7) 

m=0,l 
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where N K ^ = Fa(jp — > KKp) is the rate for kaon pair production integrated over t 
and m KK in the region of the peak. The two angular distributions Wpp and Wss 
are determined by Eq. The remaining ones are 

^5=^IpSoI«**. ^ s = -^-K°|sin^cos0, (8) 

assuming that P and S waves are produced via natural-parity t-channel exchanges, 
i.e., via Pomeron and p or u mesons, respectively ||26|| . In the general case there 
is an additional contribution to S-P interference, identical in form to the second 
term in Eq. (H) but with cos0 replaced by sin</> in Wp s and with a shifted phase 
(fig — > (fig + 7r/2 (see below). The coefficients Mjj> are given by 



\M PP \ = r^r| e -^ 



\M SS \ = r^r 2 s 



\M™ S \ = 2r^r|e-^Ree^ 

2 



< £ ±Me T ^ +iAm )5 r e -(T5-*Am)|-^ a _ i^^^-^+iAHl+^j ; (9) 



a=l 



with the upper (lower) sign corresponding to a = 1(2), and a — 1 for a = 2 and 
vice versa. The phases 4>p, m = 0, 1, arise from the elements p^o = \Pmo\ e ^ °f 
the spin density matrix. For convenience, we have also introduced the definitions 
Ar = T s — T L , Am = m L — m s , t = t 2 — h, V = V s + T L , m = m L + m s , and 
T = t x + t 2 . 

In the next section, we construct an asymmetry Ai 2 (q) to extract parameters 
sensitive to CP and CPT violation. However, we can already note here that the term 
in the production rate proportional to \M S s\ is independent of t and consequently 
is absent from the numerator of any asymmetry defined as a i-sensitive difference 
of production rates. In contrast, this term does contribute to the denominator of 
such an asymmetry, thereby making the modulation of the signal harder to see. An 
incoherent even-parity background merely renormalizes this term, further increasing 
the background but not directly influencing the signal. 

The decay rate exhibits some interesting and potentially useful properties. Con- 
sider first the case where fi = f 2 , i.e., both kaons decay to the same final state. 
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Then, it follows from the structure of the initial wave function that all terms either 
independent of or quadratic in the background (the first and third terms in Eq. (|7|)) 
are invariant under the transformation t — ► —t, whereas the linear term (proportional 
to Mps in Eq. (|7|)) changes sign. Among the consequences is that the asymmetry 
^12(9) defined in the next section contains a term linear in the background. This 
may provide a means of measuring p PS , as discussed below. If instead we consider 
the case of different decays occurring at the same time t\ = t 2 , then the terms ei- 
ther independent of or quadratic in the S'-wave amplitude are invariant under the 
transformation j\ <-> f 2 , whereas the linear term changes sign. Also, as expected, 
R(T/2,T/2,q)<x\M S8 \. 

3. Background enhancement. The parameters sensitive to CP or CPT violation are 
best extracted from singly and doubly integrated asymmetries. Define 

R(\t\,q)=R + (\t\,q)-R-(\t\,q) (10) 



and 



where 



R(q)= d\t\R(\t\,q), (11) 







^(1*1,9)= dh [ TL dt 2 5(t 1 -t 2 ±\t\)R(t 1 ,t 2 ,q). (12) 
Jo Jo 



dt x \ 
/o Jo 

Here, the time Tl is related to the acceptance for detecting Kl- Assuming a maximum 
distance of di = 5 m between the primary vertex and the detector, together with a 
boost factor 7^ ~ 5 corresponding approximately to E 1 = 10 GeV photons, we get 

TlXl = -^—? L ~ 6.5 x 10- 2 , (13) 

PkIk 

which leads to Sl = 1 — exp(— TlXl) ~ 6% acceptance for Kl detection. 
The signal can be parametrized by an asymmetry ^12(5), defined as 

a n ff L d\t\[R*(\t\,q)-R-(\t\,q)] 
U{) JdiqJ^d\t\[R+(\tlq)+R-(\t\,q)] 1 ) 

The detailed form of this asymmetry is somewhat involved. For simplicity, we adopt 
the approximations T s ~ T ~ 2Am ^> T L , which are well justified experimentally. 
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With these approximations, the coefficients in the expressions that follow hold to 
about 5%. Since for either semileptonic or 2ir decays the ratio r = 972/971 is one to 
about 0.1%, it is also useful to write r = 1 + k, where k < 1 and k can be zero for 
certain final states. Combining these approximations and definitions, we find 



ffnS L + 2 sin A0 


W PP (q) + |E^ =0 + 


Wp s (q) 


f<Pq 


(ffS L -2cosA0) 







a »® = —tj^ ^ — r A „° r , — . ( i5 ) 

Jd 2 

where A0 = </>i — 02 and 

A? = cos(0™ - 0.) + 3 sin(0™ - fa). (16) 

To gain insight into the content of the asymmetry Ai2{q), consider first decays into 
the same final state, fx = Then, r = 1, k — 0, rji — 772 = Vi an d A™ = A™ = A m . 
The expression (|T5|) therefore simplifies considerably. If the final states consist of the 
same 2n mode, then we find 

A 12 (q) r * M m "° — T . (17) 

/ d 2 q [^S L W PP (q) + ^Wss(q)} 

This provides a means to measure <fi if the background is known (or possibly a means 
to determine the background if is known). For example, the present bounds on 
CPT are limited by the precision to which <pi and 02 are known (roughly, a degree). 
Results of similar appearance arise for, say, semileptonic final states. The situation is 
also interesting for f\ = 2tt° and f'2 = tt + tt~ . Then, \r] a \ ~ 10~ 3 and r ~ 1 + 3Ree'/e. 
The issue to resolve is under what circumstances W$p 7^ enhances or masks the 
signal. 

The existing data on S-P interferometry with E 1 < 10 GeV come from two ex- 
periments, at DESY [[14]] and Daresbury [15|. Both experiments clearly identify a 



significant S"-wave background under the dominant P-wave signal. Since in diffrac- 
tion p\\ is the largest element of the spin density matrix, one would expect S-P 
interference to be dominated by the Wp s term, i.e., by p\^. However, in the limit 
t' = t — t min — > 0, the S'-wave production is dominated by nucleon helicity-flip and 
so pIq should dominate. This is consistent with the DESY results where t max ~ 0.2 
GeV. Nonetheless, the higher- 1 data from Daresbury (t max = 1.6 GeV) do indicate 
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the presence of p\°, and this despite the neglect of its contribution in the partial-wave 
analysis appropriate to the low-t behavior. Using Eq. (|T7j), it follows that a nonva- 
nishing azimuthal dependence in the S-P interference due to p\® would significantly 
enhance the sensitivity to the weak phases 0,. 

High-statistics experiments measuring all elements of the spin density matrix are 
clearly needed. An experiment of this type has recently been proposed for Jefferson 
Lab P7[ . With the planned accelerator-energy upgrade, this experiment would have 
access to approximately 10 5 7J9 — > KKp events per day (compare, for example, the 
total of 3500 events collected at DESY). Existing experimental results fail to provide 
definitive estimates for the two most relevant elements of the spin density matrix, pjg 
and p\q. The magnitude of p^ has been estimated to be of the order of a few percent, 
while no unique prediction for the phases has been found. 

The scales partially controlling the rates and the asymmetry A 12 (g) are given by 

\W PP \ ~ 0(p\\) = l, 
\Wps\/\v\ ~ Odp^DM-lWssl/lvl-Odp^M^s. (18) 

Examining the expressions given in the above expressions for the numerator and 
denominator of v4 12 (g), one discovers also a natural scale for 2ir decays given by 

a = 0(kS l T s /T l ) ~ 0(2 sin A0) ~ 0(S L \r]\T s /T L ) ~ HT 2 - 10 _1 . (19) 

The relative sizes of s and a therefore provide a separation into distinct regimes, to 
be considered in turn. 

If s < a, the asymmetry is dominated by the term 2sinA</> ~ 6Ime'/e, while 
sensitivity to the usual 3Ree'/e is suppressed due to the acceptance factor Sl- Terms 
proportional to the S*-wave background are roughly an order of magnitude smaller. 
Thus, depending on the phase of the coherent background, it may be feasible to 
extract information about <p a . 

If a < s the contribution from s dominates the numerator. As s grows, the 
denominator acquires nontrivial s dependence. This case is potentially very impor- 
tant because e'/e is no longer a factor and so the whole measurement can focus on 
improving the bounds on <p a . If an experiment can be set up in this regime, it 
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would represent a novel means of measuring two of the more elusive quantities in CP 
physics and of bounding CPT. This appears to be the most favorable case. Since 
s = O(|p 10 |)/M ~ 10 3 |p 10 |, this corresponds to a |p 10 | magnitude of greater than 
0(1CT 3 ). However, it cannot be arbitrarily large, since for s » a\ the signal be- 
comes background dominated. Although there is still information in the modulations 
via the asymmetry, the signal falls off as 1/s due to the quadratic contributions in 
the denominator, and so measurements become harder. The current experimental 
situation seems to favor values of s of the order of 0(1CT 2 ), i.e., with the K S K S 
decays from the S wave dominating the denominator. 

4- Simulation. In a real experiment, the limit of the twice-integrated asymmetry 
would not be taken, and instead the complexities of the full system would need to 
be simulated. As an initial contribution along these lines, we have implemented a 
preliminary Monte-Carlo study assuming a flux of 5 x 10 8 7/s as expected in the later 
phase of the proposed photon experiments at Jefferson Lab. This translates into 
O(10 10 ) K S K L pairs from decay per year. 

We assume for definiteness that the only nonvanishing amplitudes are the P-wave 
helicity nonflip at both photon and nucleon vertices, a™, the P-wave single helicity 
flip at photon and nucleon vertices, %, the S-wave nucleon helicity nonflip, b n , and S- 
wave single nucleon flip, \y . In terms of these amplitudes, the nonvanishing elements 
of the spin density matrix are then given by 

p\\ = pVi = 2KI 2 , Poo = 2|a'| 2 , 

p° ° = 2\bf\ 2 + 4\b n \ 2 . (20) 

We adopt the choice \a^\ ~ \b^\ ~ \b n \ = 0.1. Together with p\\ ~ 0.5, which is fixed 
by the condition 1 = trp = 2p\\ + pH + p™, this leads to Pqo ~ Poo ~ 0.01, which is 
consistent with the \ow-t DESY results. 

For simplicity, we limit attention to the case where both neutral kaons decay 
to 7r + 7r~ states. With the rates given above, the P wave contributes about 10 4 
7P — > (K°K°) p p 2(7r + 7r~)p events. Similarly, the (KsKs) s in the S wave yields 
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Figure 1: Theoretical prediction (solid) vs Monte Carlo (steps) simulations of A±2 as 
a function of the relative decay distance d as measured in the lab. 

approximately O(10 7 ) events. We adopt the former as the number of generated events. 

Under the above conditions, we have simulated R(\t\,q) and extracted the asym- 
metry A12. The result is shown in Fig. 1. The solid and histogram lines correspond 
to the theoretical prediction and the simulation, respectively. The sensitivity to the 
phase (j) a = + _ can be displayed by comparing the magnitudes of differences in the 
asymmetry simulated with <f)^ IC = 45° with theoretical predictions. In Fig. 2, these 
differences are plotted for three theoretical predictions, calculated using (fi a = <f> M 
(circles), 4> a = 4> MC + 1° (squares), and 4> a = 4> MC + 10° (triangles). Inspection of the 
two figures, in particular Fig. 2, suggests that with the anticipated number of events 
it should be possible in principle to extract the weak phases to within 0(1°) accuracy. 
With a full partial wave analysis, the sensitivity might be enhanced by another order 
of magnitude. 

Within the context of conventional quantum field theory, with CPT violation at 
the level of the standard model described by the Lorentz-violating standard-model 
extension, the effect of CPT violation on an oscillating neutral meson depends on the 
meson velocity magnitude and orientation |T7]]. However, in simulating the double- 
pion decays of the kaons, we have taken the CPT-sensitive phases to be indepen- 
dent of orientation for simplicity. This corresponds to the case where nonzero CPT- 
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Figure 2: Sensitivity of the simulated asymmetry to the weak phase <p a . 



violating phases are determined by the timelike component of the parameter Aa^ in 
the standard-model extension. For this situation, boosting the meson with a boost 
factor 7 enhances the CPT- violation effect, inducing a corresponding additive change 
in (p a by an amount 7 — 1, i.e., 4> a (l) = 70a(l)- Thus, an 0(1°) sensitivity to (fi a in a 
photoproduction experiment with boost factor 7 ~ O(10) is comparable to a O(0.1°) 
sensitivity to <p a in a similar experiment at rest. 

It would evidently be interesting to study also the prospects for enhanced CPT 
reach in the more general case of orient at ion- dependent effects. The orientation de- 
pendence leads to additional possibilities for CPT signals, including notably sidereal- 
time dependence [17 . This has recently been used by the KTeV Collaboration to 
obtain a CPT bound in the kaon system that is independent of all previous bounds 
0. In the present context, a complete analysis is likely to require more detailed input 
regarding the angular distribution of the kaon momentum spectrum and the detector 
performance. With a favorable scenario, bounds on all four independent coefficients 
for CPT violation in the kaon system could be obtained. 



5. Summary. This work has investigated some aspects of the CP and CPT sensitivity 
that could be attained by experiments involving photoproduction of neutral-kaon 
pairs. The photoproduction mechanism generates a coherent P wave as usual, but 
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also yields a coherent S wave along with an incoherent background. The magnitude of 
the coherent S wave is presently uncertain but could be significant. Under favorable 
circumstances, the resulting CPT sensitivity could be comparable to that attainable 
at conventional factories. Both analytical calculations and Monte-Carlo simulation 
indicate that interference between the P and S waves might lead to an enhancement 
of an order of magnitude in the existing CPT reach. 

We thank A. Dzierba for discussions. This work was supported in part by the 
United States Department of Energy through grants DE-FG02-91ER40661, DE-FG02- 
87ER40365 and contract DE-AC05-84ER40150, under which the Southeastern Uni- 
versities Research Association (SURA) operates the Thomas Jefferson National Lab- 
oratory (Jefferson Lab). 
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